Time Dependent Pairing Equations for Seniority One Nuclear Systems 
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When the time dependent Hartree-Fock-Bogoliubov intrinsic equations of motion are solved in the 
case of seniority one nuclear systems, the unpaired nucleon remains on the same orbital. The blocking 
effect hinders the possibility to skip from one orbital to another. This unpleasant feature is by-passed 
with a new set of pairing time dependent equations that allows the possibility that the unpaired 
nucleon changes its single-particle level. These equations generalize the time dependent Hartree- 
Fock-Bogoliubov equations of motion by including the Landau-Zener effect. The derivation of these 
new equations is presented in details. These equations are applied in the case of a superasymmetric 
fission process, that is, in order to explain the fine structure the 14 C emission from 233 Ra. A new 
version of the Woods-Saxon model extended for two-center potentials is used in this context. 

PACS numbers: 24.10.Eq Coupled channel and distorted wave model; 23.70.+j Heavy-particle decay 



I. INTRODUCTION 

In the Hartree-Fock (HF) approximation, the self- 
consistent potential for heavy nuclei is quite smooth, 
since it includes the convolution of many instanta- 
neous densities. If the potential varies in time, each 
single-particle wave function moves independently in the 
smoothly varying well. The Pauli principle is fulfilled 
being mediated permanently through the mean field po- 
tential. The two-body collisions are incorporated in the 
equations of motion only to the extend to which they 
contribute to the mean field. In principle, the time de- 
pendent HF approach treats the residual interactions ex- 
actly only if the mean field is allowed to break all sym- 
metries. A such approach leads to a very big computa- 
tional problem. In order to avoid such numerical difficul- 
ties, usually, the HF mean field is constrained to be at 
least axial symmetric. In this case, levels characterized 
by the same good quantum numbers cannot intersect. 
When the equations of motion are solved in such circum- 
stances, an individual single particle wave function will 
belong to only one orbital characterized by some good 
quantum numbers and the mechanism of level slippage 
[l| is not allowed. In general, two dynamical approaches 
are used. On one hand, the generator coordinate method 
assumes that the internal structure of the decaying sys- 
tem is equilibrated at each step of the collective motion 
||. On the second hand, the exchange between collective 
and internal degrees of freedom is neglected [|[ , so that 
adiabaticity is assumed. Usually, these behavior leads 
to the unpleasant feature that a system, even moving 
infinitely slowly, could not end up in its ground state. 
Some attempts were done to extend the time-dependent 
HF method in order to includes collision terms 0, Q. 
This problem was partially solved by introducing a resid- 
ual pairing interaction in the Time Dependent Hartree- 
Fock-Bogoliubov (TDHFB) approach || Q. This last 
method provides the possibility of level slippage for pairs 
and allows a description of the nuclear dynamics. For ex- 
ample, in the case of an even-even system, this method 
represents a tool to estimate the dissipation during disin- 



tegration processes [g, ||| . A connection with the Landau- 
Zener effect is included in the TDHFB equations. Pairs 
undergo Landau-Zener transitions on virtual levels with 
coupling strengths given by the gap A Q . Unfortunately, 
for seniority one nuclear systems the pairing residual in- 
teraction does not affect a single nucleon and, during 
the deformation of the nucleus from its initial state up 
to scission, the unpaired particle remains located on the 
same orbital. The level slippage is again forbidden for 
the blocked level. 

In the case of independent single-particles, neglecting 
residual interactions, the problem of the unpaired nu- 
cleon is solved in terms of the Landau-Zener effect. The 
Landau-Zener effect reflects a mechanism that allows the 
possibility that a single nucleon skips from one single- 
particle level to another one in some avoided level cross- 
ing regions. The probabilities that the unpaired nucleon 
arrives in different final states can be computed by solv- 
ing a system of coupled channel equations that charac- 
terizes the microscopic motion. In the following, a way to 
introduce a similar mechanism for the unpaired nucleon 
in super-fluid systems is investigated. The TDHFB equa- 
tions will be generalized in order to include the Landau- 
Zener effect. The classical TDHFB equations and the 
equations that govern the Landau-Zener effect will be 
obtained as particular cases of the new time dependent 
pairing equations. 



II. 



LANDAU-ZENER EFFECT 



The single-particle levels are function of the deforma- 
tion parameters that characterize the shape of a nucleus. 
Levels characterized by the same quantum numbers asso- 
ciated to some symmetry of the system cannot cross and 
exhibit avoided level crossings. The transition probabil- 
ity of a nucleon from one adiabatic level to another one 
is strongly enhanced in an avoided crossing region. This 
promotion mechanism is known as the Landau-Zener ef- 
feet [J. 

In Fig. [TJa) an ideal avoided crossing (J, m) between 
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III. SUPERFLUID SYSTEMS 



u m +v m a*a^ a* u m +v m a*Qs, Uj+v,a*af 




(a) (b) 



FIG. 1: Ideal avoided crossing region. Possible transition 
states in an avoided crossing region for an unpaired nucleon 
in the superfluid model. 



two adiabatic levels e k and is displayed. The diabatic 
levels are e m and tj . If the variation of the generalized 
coordinate is produced slowly and the nucleon is initially 
located on the level et , after the passage of the avoided 
crossing region, the nucleon will practically remain on the 
same adiabatic level. In this case, the motion is adiabatic 
and the nucleon follows the adiabatic state e^. If the vari- 
ation of the generalized coordinate is produced suddenly 
then the nucleon will skip with a great probability on the 
adiabatic level e k . The motion is diabatic, and the nu- 
cleon follows the diabatic state 6j. A formalism can be 
used to obtain the promotion probability. Assuming an 
n-state approximation, the wave function of the unpaired 
nucleon can be formally expanded in a basis of n diabatic 
wave functions if(r) as follows: 

#(r, t) = Ci(*Mr) exp jf e^d^j , (1) 

where the matrix elements with the diabatic states are 
£j =< tp% | H | ^ >, hij =< <pi | H | tpj >= hji, H is the 
mean field Hamiltonian and a are amplitudes. Inserting 
the wave function |T|) in the time-dependent Schrodinger 
equation, the following system of coupled equations is 
obtained: 

(2) 

Here, pi = \ci\ 2 is the probability to find the unpaired 
nucleon on the level i . 

These equations were already used to explain the reso- 
nant like structure of the inelastic cross sections in heavy- 
ion collisions [l(| E| , the fine structure in cluster emis- 
sion [T3, EH and in alpha decay H3 . the resonant struc- 
ture in the fission cross section jl5l. Il6l [r^. 



An effect analogous to the Landau-Zener one can be 
obtained by generalizing the TDHFB equations for the 
case of seniority one nuclear systems. The problem will 
be explored in the simplest possible way: a monopole 
pairing force, and a sufficiently weak pairing such that 
the nucleons are not redistributed to change significantly 
the mean field potential. In order to make the problem 
tractable, two approaches are investigated. The first one 
is valid for a low lying levels system with a small number 
of avoided crossing regions, so that variations of densities 
Pi and pairing moment components ft; due to the blocked 
level can be neglected. The second one takes into account 
the blocking effect, that is, the fact that pn^n) and Ki( m ) 
depend on the blocked level m. 



A. Low lying levels 

Using quasiparticle creation and annihilation operators 
a~ji and aj. 

a k = u k a k -v k a^] 
a k = u k a k + v k a k ; 

a k = u k<4 - v* k a k , (3) 
a k = u k a^ + v* k a k ; 

it is possible to construct some interactions that help us 
to promote the nucleon from one diabatic level to an- 
other. The two situations plotted in Fig. Q]can be mod- 
eled. In the plot[lJa), the single particle follows the dia- 
batic level ej while in QJb) it stays on the adiabatic one 
et. Here a k and a k denote operators for creating and de- 
stroying a particle in the state k, respectively. The state 
characterized by a bar signifies the time-reversed partner 
of a pair. The parameters v k and u k are the occupation 
and vacancy amplitudes, respectively. Because only the 
relative phase between the parameters u k and v k mat- 
ters, in the following u k is considered as a real quantity 
and v k a complex one. The interaction able to promote 
the unpaired nucleon from one adiabatic level to another 
must be given by product of operators of the type © . 

In order to obtain the equations of motion, we shall 
start from the variational principle taking the Lagrangian 
as 

d 

5L = 5 < w\ H -ih— + H' - XN I ip >, (4) 
at 

and assuming the many-body state formally expanded as 
a superposition of n time dependent BCS seniority one 
diabatic wave functions 

n 

<p(t) = Y, c>n(t)a+ H («,(t) + v l (t)a+a+) . (5) 
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The Lagrangian contains several terms. The first one is 
the many body Hamiltonian with pairing residual inter- 
actions 

H(t) = e k (t)(a+a k + a^a k ) - G a+a^aiaj. (6) 

fc>0 k,l>0 

The residual interactions between diabatic levels charac- 
terized by the same quantum numbers that are responsi- 
ble for the Landau-Zener effect are assumed on the form: 

= E"jy« hij (t)(uidf - v?Oi)(ujaj - Vjat). 

The sum runs over diabatic levels i and j. The number 
of particle operator is: 

n = + a t a ~k) ( 8 ) 

fe>0 

After some calculations, as detailed in Appendix \X[ the 
next time dependent coupled channel equations are ob- 
tained pi]: 

n 

ihpi = ^2p m {niA m - «*A m } , (9) 

m 

n 

ihki = {(2pi - 1) A m + 2ki (e, - A)} , (10) 

m 



ifvPm — / y ^mj {Smj $jm) i (H) 



the operation is realized on the n possible diabatic states 
of the unpaired nucleon. In this paper, the sum over 
pairs energy generally runs within the index k. When the 
single-particle sum over k is realized only for one part- 
ner of each reversed pair the result is multiplied with a 
factor 2. The index k runs over a workspace that allows 
the pairing force to operate only over a a finite number 
of active levels around the Fermi energy. 

B. Blocking effect 

If the blocking effect is taken into consideration, each 
seniority one BCS wave function is characterized by its 
own set of p and k values, and the trial wave function is: 

n 

<p(t) =^2c m (t)a+ Yl (ui(m)(t) + vi( m )(t)afaf) . 

m l=jim 

The Landau-Zener interaction is postulated as follows: 
xa t(j) a m llfe^ij a k(j)at <*ka+ {i) 

(15) 

where the quasiparticle creation and annihilation opera- 
tors 

a k {j) = u kU) a k - v kU) ap 

a k(j) = u k(j)a-k + v k( 3 )a k ; 

a t{j) = u m a k ~ v k(j) a k' (16) 

a t( 3 ) = u k{j) a t + v l{j) a ^ 



ihS jm = S jm {-£ (| A m | 2 - I Aj \ 2 ) + (e m (t) - £j (t)) 



2K m + ^-2Kj)j2?Pl^l} 



+ Td^m,j [hml(t)Sjl - hjl(t)Sl m ] + h m j(t)(pj - p m ) 

The next notations are used: 



(12) 



A, 

a: 



= GJ2 
= GE k 



k^m 



k^m "-fci 



Kk; 



K k 

Pk 

Pm 
S 



UkVk; 

\v k | 2 
\c m | 2 



(13) 



p are single particle densities, n are pairing moment com- 
ponents, and p m denotes the probability to have an un- 
paired nucleon on the level m. p and p are real quan- 
tities while k and 5* are complex ones. In analogy with 
k, S (having the property | Sj m \ 2 =PjP m ) can be called 
as unpairing moment component. Whenever the upper 
limit n is specified for a sum, it is implicitly assumed that 



are now associated to each blocked level (j). This inter- 
action considers only an approximate way to describe the 
full phenomenon. In this context, if a diabatic wave func- 
tion i is "reflected" in an avoided crossing region 
this wave function is transformed in a component of the 
diabatic wave function j. The reality is more compli- 
cated. When a diabatic wave function i is "reflected" in 
an avoided crossing region, this wave function must be 
split into two parts: a transmitted diabatic wave function 
i and a reflected adiabatic wave function j' . That means, 
the number of wave functions must be doubled after the 
passage of each avoided crossing region. Therefore, treat- 
ing the more realistic situations, the system of coupled 
channel equations becomes much more complicated. For 
simplicity, in our approximations we considered only a 
superposition of n diabatic wave functions, that means, 
the diabatic wave function i is forced to contribute to 
the amplitude of the diabatic wave function j (which is 
not always equivalent to j'). After some calculations, as 
detailed in Appendix [B] a new set of TDHFB equations 
results: 



iftpl(m) = «i(m)A£ 



K ;*(m)Am, 



(17) 
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iftki(m) = i?Pl(m) - l) A m + 2K ;(m) (ei - A m ) , (18) 



where the solutions of the TDHFB equations are used. 
In the frame of our model, the difference 



(19) 



ifo$jm — Sj m ^ G (l 



IE 



+(e, 

(a, 



>(£)-£.,(£) -A r , 



v fe(m) 



- I A, | 2 ) 

2 



Efc 



A J K fc(j) 



+ EL,j [h m l{t)Sjl 



hjl(t)Sl m 



Pfc(j) 



- 1 



h mj(t)(pj 



1 



(20) 

The same notations as in the previous approach are used. 

Two main differences arise between Eqs. (jiJjl- CH?]) and 
(fT7| - (|20p that are implicitly determined by the hypothe- 
sis assumed in their derivation. Firstly, in Eqs. ([9|)- (fr2|) 
the values of p and k are obtained through a weighted 
sum that runs over unpaired states while in Eqs. (JT7J)- 
(|T8]) these quantities belong to only one diabatic wave 
function. Secondly, in Rel. (|20[) . Sj m depends on all 
the densities p and pairing moment components k of 
the implied two diabatic wave functions j and m. As 
a consequence of these differences, the number of dif- 
ferential equations increases n time in Eqs. (fT7|) -(fT8 |) . 
Another consequence is that 2j2 k _ £m Pk( m ) — N — 1 for 
each diabatic wave function m in Eqs. (fT7l) - (fT8| while 
2 Efc Pk =N + 2p F -l in Eqs. [j9"])-([T2"]) where p F denotes 
the single-particle density of the Fermi level in the initial 
ground-state configuration. Finally, the chemical poten- 
tial A has values associated to the diabatic state under 
consideration in Eqs. (fl"7]) -(fl"8 ]) . 



IV. ENERGY 

In this section, only the equations associated to the 
blocking level approach are displayed. For the low lying 
level approach, the index (m) must be dropped. The 
ground state energy Eq of any deformation is obtained 
in the framework of the BCS formalism by considering 
the Fermi level ep populated with the unpaired nucleon: 



k^F 



E F -G 



k=£F 



G X! P k(F)- 

k^F 



in the static, lower energy state. For the same deforma- 
tion, the energy of an adiabatic state m is obtained by 
considering the unpaired nucleon located on the diabatic 
state under consideration: 



E m — 2 ^2 Pk(r, 
k^rn 



,-G 



Y K k{ m ) 
k^m 



k^rn 

(22) 



A£L 



E„ 



En, 



(23) 



behaves as a specialization energy. So, as inferred in Ref. 
[Ill the quantity AE m must increase the potential bar- 
rier tunneled by the nuclear system. Different barriers 
are obtained for each diabatic state under consideration. 
These appear as dynamic excitations during the decaying 
process. Combining excitations with occupation proba- 
bilities of diabatic states, we obtain 



(24) 



for the average energy and 



AE = J2Pm&E m , 



(25) 



for the averaged dissipated energy during the decay. As 
mentioned in Ref. [201 ] . the collective kinetic energy is 
temporarily stored as a conservative potential. This en- 
ergy subsequently decays partially to the dissipation. 

The equations @-||T2J) and (|T7 1) -([2"0" ]) involves only 
single-particle energies. They conserve the average num- 
ber of particles because 2 J^k^m Pk(m) — N — 1 for any 

TO ( or 2 Y,kPk = N + 2 Pf - 1) and X^Pm = 1- The 
average energy can evolves in time as follows: 



E - El Pm {2 J2k^m Pk( m )£k + e m 

@ I Efc/m K k(m) \ 2 —GJ2k^m Pk(m) } • 



(26) 



For a stationary system, for which e=0 and G—0, the to- 
tal energy is conserved, even if individual values of p, p 
and k may still be varying with time. In our treatment, 
the chemical potential has the values A m obtained from 
BCS equations for each energy levels workspace associ- 
ated to the diabatic wave function m. 



V. GENERALIZATION 

If the blocked levels are eliminated, the system (f9"])- (Ti"2")) 

reduces to: 

ihpi = K t A* ~ k\A, , 
ihki = (2 Pl - 1)A - 2k, [ei(t) - A(t)] , 

the well known TDHFB equations @, Q. On another 
hand if the pairing is neglected, the third equation of the 
system (]TT]) can be written 



ih (c rn C m -f- C m C m ) — ^ ^ hmj (pjC m -{- CjCjj 



(28) 
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Introducing explicitly the time dependence of the ampli- 
tudes c m 



c m (t) = c 0m {t)exp 



the next relation is obtained: 



e m (r)dT , (29) 



c 0j c 0m exp 



CojCom exp J *( ej - e m )dr 



(30) 



The last relation is an equivalent form of the Landau- 
Zener Eq. ([3]) obtained in the frame of the single-particle 
model. Furthermore, if the pairing interaction is ne- 
glected, p and k can be either zero or one, and the fourth 
equation of the system ([12]) reduces to 



+ Si^m.j (hmlC*Ci — hjiC^C m ) + hj m (c*Cj — C^Cm) ■ 

(31) 

After introducing the exponential dependence, the next 
relation emerge: 



ih (c 0m c^ + CojCom) exp (-jr J (e m - ej)dr) 
= Ej^mj [ h mlc* oj Coi exp (-i J(e t - e^dr) 
- hjic^comexp (-7- / (e m - ei)dr)] 

+h m j (cg^Coj — c5 m CQ m ) . 



(32) 



That is another form of the Landau-Zener relation ([2]). 
So, the Landau-Zener equation for single particle systems 
(without residual interactions) and the TDHFB equa- 
tions for quasiparticles are two particular cases of the 
coupled channel equations (|9"))- (fT2"j) . So, this system rep- 
resents a generalization of the TDHFB equations in the 
case of seniority one nuclear systems. Similar arguments 
are valid also for the system (Tr7|) - (f2"0"|) . 



VI. RESULTS 

To solve the TDHFB equations, only the variations of 
the single-particle energies et are needed. The simplest 
way to obtain the evolutions of single-particle energies is 
to consider a time-dependent single particle potential in 
which the nuclcons move independently. As evidenced in 
Ref. @ , such a description is within the spirit of the more 
rigorous Hartree-Fock approximation, which defines the 
potential self-consistently. 

The 14 C emission from 223 Ra will be treated. The 
fragments issued in this reaction are spherical while the 
parent is little deformed, allowing a description in terms 
of a nuclear shape parametrization given by two spheres 
smoothly joined within a third surface. 

A fine structure in the 14 radioactivity of the 223 Ra 
was observed in 1989 [U, [22|, [23j]. In the first exper- 
iment, the results indicates that 15±3 % of 14 C decays 
are transitions on the ground state of the daughter, while 



81±6% are transitions on the first excited state. In Ref. 
[24| . using the M3Y potential, it was evidenced that the 
preformation probability must be more favorable for the 
excited state than for the ground state with a factor of 
180. Such a value cannot be accounted from theoretical 
models [25[ without taking into account dynamical in- 
gredients. This is the main reason that the fine structure 
phenomenon was selected to validate our equations. 

The deformation energy of the nuclear system is the 
sum between the liquid drop energy and the shell ef- 
fects, including pairing corrections. The macroscopic en- 
ergy is obtained in the framework of the Yukawa-plus- 
exponential model extended for binary systems with dif- 
ferent charge densities [26| . The Strutinsky prescriptions 
[27l | were computed on the basis of a new version of 
the superasymmctric two-center shell model. This ver- 
sion solves a Woods-Saxon potential in terms of the two- 
center prescriptions as detailed in Appendix [Cl 

Because the pairing equations diverges for an infinite 
number of active levels, a limited number of levels are 
used in the calculations: 31 levels above and 31 levels 
under the the unpaired Fermi level in the initial ground 
state configuration that is, N — 1=62. These levels are 
selected in terms of the spin projection Q on the symme- 
try axis and kept as a single particle energies workspace. 
A constant value of the pairing parameter G — 0.13 MeV 
is used. 

The least action trajectory was obtained by general- 
izing in a three-dimensional space the method initiated 
in Ref. 12811 a nd used extensively to describe the fission 
processes [l5l Fill \vj\ . The inertia is computed within 
the Werner- Wheeler method. The trajectory of the de- 
caying system is obtained simultaneously as function of 
three generalized coordinates, that is, the elongation R 
(the distance between the centers of the nascent frag- 
ments), the necking parameter C = S/R3 (the curvature 
of the intermediate surface) and R1/R2 (the ratio be- 
tween the radii of the heavy fragment R\ and that of 
the light one i?2)- These parameters are explained in 
Appendix [Cl In Fig. [2Ja) the deformation energy V 
of the nucleus is plotted as function of the elongation 
R. Three excitations of the nuclear systems that corre- 
spond to three adiabatic wave functions are also plotted 
with dotted lines. These excitations are added to the 
deformation energies obtained in the framework of the 
macroscopic-microscopic model in order to calculate the 
penetrabilities as show below. In Fig. [U (b) and (c), the 
variations of the necking and mass-asymmetry general- 
ized coordinates are displayed. At R ss 10 fm, a system 
formed by two spherical tangent nuclei is obtained. The 
Woods-Saxon potential is presented in Fig. [3] for a se- 
quence of nuclear shapes along the least action path. 

The single-particle level schemes for neutrons and pro- 
tons must be computed along the minimal action trajec- 
tory in order to solve the time dependent pairing equa- 
tions. It is known that 223 Ra has the spin | emerging 
from lin/2- Adiabatically, the unpaired neutron reaches 
the 2gg/ 2 level of the daughter 209 Pb. As also evidenced 
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FIG. 2: (a) Deformation energy V as function of the distance 
between the centers of the nascent fragments R. The three 
excitations due to the diabatic levels £j, i — 1,3 are also 
plotted with dotted lines, (b) Variation of the curvature of 
the median surface and (c) of the mass-asymmetry parameter 
as function of R. 



FIG. 4: Neutron energy diagram along the minimal action 
path as function of the distance between the centers of the 
fragments R. The levels with spin projection S7 of interest 
are plotted with thick lines. The levels are labeled within 
the spectroscopic factors. At the right the first column corre- 
sponds to the daughter nucleus while the second one is related 
to the 14 C. 
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FIG. 3: Mean field Woods-Saxon potential Vb as function of 
the cylindrical coordinates p and z for different values of the 
elongation along the minimal action trajectory, (a) Elongation 
R=2 fm; (b)i?=5 fm; (c) R=W fm; and (d) R=15 fm. 



FIG. 5: (a) Selected neutron energy levels that can be occu- 
pied by a single neutron as function of the internuclear dis- 
tance. Thick lines are the diabatic levels e», i=l-3, while thin 
lines are used for the adiabatic ones, (b) Interactions energies 
hij in the avoided crossing regions. 
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FIG. 6: (a) Occupation probability pi for the diabatic level 
ei as function of the internuclear distance R. (b) Occupation 
probability p2 of £2. (c) Occupation probability p$ of £3. (c) 
The solid line corresponds to the dissipated energy AEi for 
ei, the dashed line is the dissipated energy AE2 for £2, the 
dot-dashed line is the dissipation energy AE3 for £3, and the 
dotted line is the dissipation energy AE P for the proton level 
scheme. 



Our goal is to compute the occupation probabilities of 
the 3 levels of interest at the end of the disintegration 
process. For this purpose, Eqs. (fTS|) - (f2T))) are used. 
Some features concerning the less rigorous low lying lev- 
els approach ([9|)- (fT0)l can be found in Ref. [l8| 

In Fig. O the three selected diabatic levels e m (m = 
1, 2, 3) are plotted together with the interaction energies 
hij determined by using spline interpolations around level 
crossings. Diabatically, the unpaired neutron, initially lo- 
cated on the level ei that starts from the spherical orbital 
1*11/2 will arrive on the final state lin/2, that is the first 
single particle excited state of the daughter, after the 
passage of three avoided level crossing regions. 

The initial conditions are determined by solving the 
BCS equations for the 3 possible seniority one wave func- 
tions at R sal. 5 fm, where the first minimum of the de- 
formation energy is located. The time dependent pairing 
equations are integrated numerically using the Runge- 
Kutta method. The occupation probabilities p m and the 
dissipated energies given by formula (|23p arc determined 
along the minimal action path for a internuclear velocity 
^ = 1.4 x 10 6 m/s. This value can be translated in a 
time required to penetrate the barrier of about 1.4 x 10~ 20 
s. 

In Fig. [S]the probability of occupations p m with an un- 
paired neutron of the three diabatic levels are presented 
as function of the internuclear distance. In the bottom 
panel, the three dissipated energies AE m and the energy 
dissipated in the proton level scheme arc also displayed. 

The branching ratio r between the partial half-life for 
transitions to the ground state of the daughter and the 
partial half-live to the first excited state is given by 



pi exp(-Ki) 
p 2 exp(-if 2 )' 



(33) 



where the index corresponds to the diabatic level e\ or 
62, and 



K m = - 



\J fi[V(R) + D m (R) - V(R ss )]dR (34) 



in Refs. [HI, [T^|, the fine structure in the 14 C radioactiv- 
ity can be understood by an enhanced transition prob- 
ability of the unpaired neutron from the adiabatic level 
51 = 3/2 emerging from lin/2 to the adiabatic level with 
the same spin projection f2 that emerges from lj'15/2, m 
terms of the Landau-Zener effect. The level scheme of 
Fig. Q]shows that (adiabatically) for fl = 3/2 the l«n/2 
level reaches the 2g g / 2 daughter state, the lj'15/2 level ar- 
rives on the lin/2 one - I n this respect, the level scheme 
calculated within the Woods-Saxon model is in qualita- 
tive agreement with that obtained within the modified 
oscillator model [13]. The O = 3/2 levels subjected to 
avoided level crossings that can give rise to Landau-Zener 
effect are plotted with thick lines in Fig. |4j Two adjacent 
levels with 57 = 3/2 are also plotted with point-dotted 
lines to show that no other avoided level crossings are 
possible if the unpaired neutron originates from Xin/2- 



are the WKB integrals. Here i? gs is the ground state 
elongation, R m is the exit point from the barrier for the 
channel m, V(R) is the macroscopic-microscopic energy, 
D m (R) = AE m (R) + AE p (R) is the dissipated energy, 
AE m being the specialization energy given by Rel. (|23[) . 
AE p denoting the dissipated energy of the proton subsys- 
tem, and /i is the reduced mass. AE p is calculated with 



Eqs. [[27), The barriers obtained for V(R) + D m (R) are 
plotted in Fig. [2] (a) with dotted lines. The experimental 
values of r range between 5.4 and 5.9. Our theoretical 
value is r=5, which is in an excellent agreement with 
experimental data. 

In conclusion, two approaches that generalize the 
Landau-Zener equations for seniority one superfluid sys- 
tems are presented. The new formalism is valid for 
any kind of mean field approximations that include a 
monopole pairing field. The equations that describe 
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our approaches offer information about the spectroscopic 
amplitudes and the dissipated energies in different final 
channels. The new equations were used to reproduce the 
qualitative and quantitative features of the fine struc- 
ture phenomenon in cluster decay. Up to now, this phe- 
nomenon was not described adequately in the frame the- 
ories that do not include dynamical ingredients, as ev- 
idenced in Ref. [29| . Within the time dependent pair- 
ing equations, a good agreement with experimental data 
was obtained. A new version of the superasymmetric two 
center shell model based on a Woods-Saxon potential was 
developed and used in this context. 
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APPENDIX A 



where | (f> k >— | >, due to the normalization. The 
high order term | v k | 4 of formula (|A3| is neglected. The 
equality ii k = — {v k v k + v k v k ) / (2u k ) is also used. 

To minimize the functional, the expression (|AT|) is de- 
rived with respect to the independent variables and 
vi. Two equations follow: 

Em I c m | 2 {2v*(e t - A) - G[£ fc/m n k (-^f ) 



Em I ^ | 2 {2v l (e l - A) - G[£ Mm nl (-^ 



0. 



(A7) 



(A8) 



where the notations for densities p =| v | 2 and pairing 
moment components k — uv are introduced. 

The condition of conservation of the number of parti- 
cles 



Vk 



N + 2 PF - f, 



(A9) 



The TDHFB equations when the blocking level is ne- 
glected will be derived in this Appendix. Following the 
same prescriptions as in Ref. @, using the time depen- 
dent Hamiltonian © within corrections ([7]) , and the trial 
wave functions (O, the expected value of the Lagrange 
function is obtained: 



was used so that 



= E r ' 



<ip\H - ih-^j + H' - XN | tp > 



{2E Mm K I 2 (e* - A) + (e m - A) 

~ G I Efc^m U t V k ?}- 
lfi Em I C m I 2 iYtk^rn U V kVk ~ V* k V k )} - iHYsZ^rnCn 



m,j=jtm hmjC m c j 



The following identities were used: 



(Al) 



< <P I J2k e k(a+a k + a±a k ) \ ip > 
= Em I Cm | 2 [e m + 2 J2 k ^ m £ fc I Vk I 2 ]; 



< V I Efei <4 a t a l a l I V > 

Em I Cm | 2 (Efc^m I v k | 4 + Ei^m «J«* Efe#m M fc w fc 

(A3) 



< | afaj + a+cti | <p >= c*Cj + c*c. t ; 



and 



= E' 



<v\m\ ( P> 



C* m C m + | Cm | 2 Ek/ m («fc«l= + V * k Vk) 

because, as evidenced in Ref [9(, 

bUk + <P*Md 3 r = 0, 



(A4) 



(A5) 



(A6) 



^(v* k v k + v* k i k ) = 0, 



and 



1_9 

2du 



(A10) 



(AH) 



(A12) 



Multiplying Eqs. (TAT)) and (|A"8)) with u ; and u*, re- 
spectively, and subtracting, the first TDHFB equation 
© is obtained: 



iHpi = 



where J2 r 



Em I Cm | 2 i^A m -<A ro } ^ 
Em I Cm | 2 

1 and A m = GJ2 k ^ m K k- 



(A13) 



Another equation can be obtained: 

^--^ + ^ = -^r^E:icml 2 

x{ Ki A^- K ;A m }+^Em|Cm| 2 

\AU) 

so that the second TDHFB equation (fT0|) follows: 
ififtj = J2 I Cm I 2 {(2p; - 1) Am + 2m (ei - A)} . (A15) 

m 

Using the property 

n 

£ I c™ | 2 = 1; (A16) 
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so that 



(A17) 



the Eq. (|Aip is derived with respect c m and c* m and set 
to zero. The next relations are obtained: 

-ihc* m = c* m {2 £ fc#m I Vk | 2 (e k - A) + (e m - A) 

~G | Ek^n U >k | 2 ] 

(A18) 

ihc m = C m [2 Efc^m \ V k? (e*! ~ A) + (e m - A) 

-G|E M ™<^I 2 ] „ (A19) 

- ific mE^ m |( U k** ~ "fe V *)] + Ej^m hmjCj, 



Multiplying the relations (|A18| and (|A19|) with c m and 
c* m and subtracting them the next relation follows 

n 

ih(c m c m + c m c m ) = h mj (c 3 c* m - c*Cm). (A20) 

jjtrn 

It is a form of the third TDHFB equation (jTTJ) . For a 
passage through only one avoided crossing region (to, j), 
only two amplitudes, c m and Cj can change. On another 
hand, from the condition of conservation it can be ob- 
tained: 



-(c iC *+c,-c*) (A21) 



This condition is fulfilled by the above Eq. (|A20|) . so the 
equation conserves the norm. Changing indexes, multi- 
plying with amplitudes and subtracting relations (|A18[) 
and (|A19[) the next equation follows: 

ih{t m c* + c*c m ) = c m c* | A m | 2 - | A, | 2 )+ 

2 | Vj | 2 ( £j — A) + (e m — A) - A m 
-2\v m | 2 (e m - A) - (ej - A) - Aj} 

+C m C* l -^{v* m V* m - V* m V m - Vpj + V*Vj) 

(A22) 

From Eqs. (|A7|) and (IA8|) . the expressions in the last 
parenthesis that involves v and v* can be obtained: 



f(vfv l -vfv l ) = 2p l (e l -X) 



(S-«r) + s^ (#-«■) 



(A23) 



Using the notations (|13[) and rearranging the terms the 
fourth TDHF equation (fT2")) is obtained. 



APPENDIX B 

The TDHFB equations when the blocking effect is 
taken into consideration are derived in this Appendix. 



Using the corrections (I15p . and the trial wave functions 
(fT3| . the expected value of the Lagrange function is: 

<ip\H- ih§- t +H'-XN\ip> 

= Em I C m I 2 { 2 Efc#m I V k(rn) ? Ofc ~ Am) 
+ (fm - A m ) - G | Efc^m u fc(m)«fc(m) | 2 } 
-^Em I Cm HEfc/m IK^^M ~ '"l( m ) V k(m))] 
~ *^Em C m c m + Era.j/m h m jC m Cj 

(Bl) 

In order to minimize the functional, the expression (|F31[) 



is derived with respect to the independent variables v 



l{m) 



and f;( m ) by taking into account the subsidiary condition 
(|A12|) , The next relations follows: 



-G 



J Cm | 2 {^*i {m) {ei - X m ) 

V Kur > ( V " {m)V " { 



fc/m "'fc(ro) 



(m) 

) 



l{m) 



= 0, 



(B2) 



Em I C m I {2f/( m ) (e; - A m ) 



-G 



Ul(m) 



E 



fc^m K fe(m) 



/ ^(m)^i(m) \ 

V 2 «k™) J 



Efc^m K k(m) 



ihv, 



l{m) 



= 0. 



(B3) 

This system can be solved by considering that the ex- 
pression in the curly bracket is zero for each value of to. 
Following a similar way as in Appendix 1, the first two 
TDHFB equations associated to an unpaired nucleon in 
the state to emerge: 



ikpi( m ) = «;(m)A m - K,* (m) A m , 



(B4) 



ihki( m) = (2p /(m) - l) A m + 2K; (m ) (e/ - A m ) ■ (B5) 

To obtain the probability that an unpaired nucleon is 
located on a state to, the expression (|B1[) must be derived 
with respect c m and c* m . Two equations follow: 



-ihct 



2 Efc^m I V Hm) I 2 (ffe ~ A m ) + E m - A, 
^ G I Efc#m U fe(m)Wfe(m) | 2 



fe#m 2 V IJ fe(m) W fe(™ 



k(m 



)«fe(m)) 



+ Ei^rn hrnjC* — 0, 



(B6) 



c m 2 E^fc^m I w fc(m) | 2 ( e fe — A m ) + £ m ' 
~ G I EMm U fe(m) U fe(m) P 
" ihc m J2 k ^m lK( m )*fc(m) ~ «fc(m) W fc(™)) 
+ Ej/m ^rojCj = 0. 

(B7) 

Multiplying with complex conjugates and subtracting, 
the next relation 



En 



( C 3< 







Cm) ) 



(B8) 
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is obtained. Using the notations (| 13[) the Eq. ([19]) fol- 
lows. 

From relations (|B6j) and (|B7|1 , another relation can be 
deduced 

ih(c*c m + c m c*) 
= c m c* [-i (I A m | 2 - | | 2 ) + (e m - ej - A m + Ay) 

+ 2 Efc^m Pfc(m) (Cfe - A m ) - 2 p fc y) (ej, - Xj 

— Y c mS L2-»A^mk u fc(fn) U *("») _ V k(m) V H™)) 

- Efe^K (i )*fcO-) -*fc(j)^(j)) 
+ Si^m hi m cic* — J2l^j hljCfCm. 

(B9) 

The derivatives £j and i>* appear in the previous expres- 
sion. In order to evaluate quantities where these deriva- 
tives intervene, Eqs. (|B2|) and (|B3[) are used. The next 
relation follows. 

i h 



— {(2«i( m )(ei - A m ) 



-rrv k* ( °'i-)'"' m| ''"'"' 

+ ( U !H - £fc#m «fc(m)]} 

+^{(2^ (m) (Q-A m ) 



J2k^rn K *k(m)\} 



Pl(m)( e l - <^m) + 



+Pl(m){ e l - ^ 



pit" 

2 2k* . ^/(m) 2K l{m) > 2 



2 2K !(m) 



2p;(m)(Q — A m ) 



Kj(m) I + — I 7^-; ~ K l( m ) 



(BIO) 



by the rotation of a circle around the axis of symmetry 
as displayed in Fig. [7] The parametrization is character- 
ized by 5 degrees of freedom that can be associated, for 
example, to the elongation (R — Z2 — z\), to the necking 
(C = S/R3), to the mass asymmetry (rj = 01/02), to 
the deformations of the two fragments (pi/di, i = 1,2). 
Treating the 14 C emission, the deformations of the two 
fragments can be neglected and the mass asymmetry pa- 
rameter is considered as i] = i?i/i?2- The mean field po- 
tential is defined in the frame of the Woods-Saxon model: 



Vo{p, *) = — 



1 + exp 



A(p, 2 ) 



(CI) 



so that relation (IB9I) becomes 



where A(p, z) represents the distance between a point 
(p, z) and the nuclear surface. This distance is measured 
only along the normal direction on the surface and it is 
negative if the point (p, z) is located in the interior of 
the nucleus. V c is the depth of the potential while a 
is the diffuseness parameter. In our work, the depth is 
V c = Vqc[1 ± k(N - Z )/N + Z )] with plus sign for pro- 
tons and minus sign for neutrons, Vq c — 51 MeV, a=0.67 
fm, k=0.67. Here A , N and Z represent the mass 
number, the neutron number and the charge number of 
the parent, respectively. This parametrization, referred 
as the Blomqvist-Walhlborn one in Ref . , is adopted 
because it provides the same radius constant ro for the 
mean field and the pairing field. That ensures a consis- 
tency of the shapes of the two fields at hyperdeforma- 
tions, i.e., two tangent ellipsoids. 

In Fig. [3] the mean field potential Vq is plotted as 
function of cylindrical coordinates p and z for four nu- 
clear shape configurations obtained along the minimal 
action path. 

The spin-orbit coupling is assumed of the form 



ifcSjm — C m Cj- { G (I Am I 

+ ( e m — £j — A m + Aj) 



' 2 X/fe^m 
+ 2 Efc#j 



A, 



A, I 2 ) 



A,- I SiSH - /,- 

J 1 K k(j) 



Ho) 



a: 1 ^-K fca) 



(Bll) 

After some rearrangements of terms, and using notations 
([15]) the Eq. (2D]) is obtained. 



APPENDIX C 



-2A 



1 



2mc 



(W x p)s 



(C2) 



where A=35 is a dimensionless coupling constant, m is 
the nucleon mass while c denotes the speed of the light. 
The spherical components of the operator 



L = W x p 
in cylindrical coordinates are 



L ± = The ± lv (Wo ± _ dVo d_ ± .mi d_ 
\ dp dz dz dp dz p dip 



(C3) 



(C4) 



A two-center shell model with a Woods-Saxon poten- 
tial was developed recently [31| . An axial symmetric 
nuclear shape parametrization is used to determine the 
mean field potential. This nuclear shape parametrization 
is given by two ellipsoids (of different semi-axis and ec- 
centricities) smoothly joined with a third surface given 



so that 



Ls 



dV l d 
dp p dip 



■ L z s, 



(C5) 



(C6) 
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The next step is to obtain the solutions of the 
Schrodinger equation 

-^A + V (p, z) + V ls (p, z) + V c {p, z)] z, <p) 

= E9(p,z,<p) 

(C7) 

For protons, a Coulomb term Vc is added as in Ref. 
[30] | . No analytical solutions can be found for such poten- 
tials. A suitable eigenvector basis able to diagonalize the 
Woods-Saxon potential can be obtained with the double 
center harmonic oscillator model. 

A complete analytical eigenvector basis can be only ob- 
tained for the semi-symmetric two-center oscillator. This 
potential corresponds to a shape parametrization given 
by two ellipsoids that possess the same semi-axis perpen- 
dicular on the axis of symmetry. The potential is 



V (p,z) 



- J 2 



\m<J 1 zl {z — ci) 2 + irru' 



\muj 2 z2 (z 



C2 



1 2 
— TtlLU 

2 p ' 



z < 0, 
z > 0, 



(C8) 



where uj denotes the stiffness of the potential along dif- 



ferent directions as follows, 



w 



= 41A 



-1/3 



A 1 / 3 



UJ - 



R<\ 



a 2 ' 

in order to 



ensure a constant value of the potential on the surface. 
The origin on the z-axis is considered the location of the 
plane of intersection between the two ellipsoids. 

An analytic system of eigenvectors can be obtained for 
Vq by solving the Schrodinger equation: 



2mg 



A + V (p,z) 



*(p,z,<p) = E*(p,z,<p) (C9) 




The analytic solution of Eq. 

ansatz 



C9P is obtained using the 



with 



9(p,z,<p) = Z(z)R(p)*(<p) 



$ m (<p) = — = exp(zm^)) 

V Z7T 



(CIO) 



(Cll) 



FIG. 7: Nuclear shape parametrization. Two intersected el- 
lipsoids of different eccentricities are smoothly joined with a 
third surface. Two cases can be obtained: (a) the curvature 
of the circle of radius R3 is positive (s=l) and (b) the curva- 
ture of R3 is negative (s=-l). The elongation is given by the 
distance between the centers of the ellipsoids R = Z2 — zi 



Rnmip) = \lJ^my apeXP \~^2~) ^ mL n( a > 2 ) 

(C12) 

f C Vx exp(- a ^ ( 7 Cl)2 )H yi [-a zl (z + ci)], 

Z (z)=\ Z < °' 

I C V2 exp (- <^ z - c " f ) Mz - c 2 )}, 

i z > 0, 

(C13) 

where L™(x) is the Laguerre polynomial, H„(() is the 
Hermite function, 014 — (rriQUJi/h) 1 ^ 2 (i = zl,z2,p) are 
length parameters, and C Vi denote the normalization 
constants. The quantum numbers n and m are integers 
while the quantum number v along the z-axis is real and 



has different values for the intervals (— 00, 0] and [0, 00). 
Imposing conditions for the continuity of the wave func- 
tion and its derivative, together with those for the sta- 
tionary energy and orthonormality, the values of V\, ^2, 
Cjy^ and C1/2 

are obtained. Details concerning these solu- 
tions and expressions for the normalization constants are 
found in Refs. [HI, [Hj]. For reflection-symmetric shapes, 
the solutions along the z-axis are also characterized by 
the parity as a good quantum number. The basis (|C13[) 
for the two-center oscillators can be used for a various 
ranges of models more of less phenomenological [3J, |35| . 
On another hand, they are different ways to obtain the 
single-particle energies for a two-center Woods-Saxon po- 
tential. Other recipes are given in Ref. (36[ where the 
potentials are expanded in terms of harmonic oscillators 
functions. 
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